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Abstract 

To enable arriving aircraft to fly optimized de- 
scents computed by the flight management system 
(FMS) in congested airspace, ground automation must 
accurately predict descent trajectories. To support de- 
velopment of the predictor and its uncertainty models, 
descents from cruise to the meter fix were executed us- 
ing vertical navigation in a B737-700 simulator and a 
B777-200 simulator, both with commercial FMSs. For 
both aircraft types, the FMS computed the intended de- 
scent path for a specified speed profile assuming idle 
thrust after top of descent (TOD), and then it controlled 
the avionics without human intervention. The test ma- 
trix varied aircraft weight, descent speed, and wind 
conditions. The first analysis in this paper determined 
the effect of the test matrix parameters on the FMS 
computation of TOD location, and it compared the re- 
sults to those for the current ground predictor in the 
Efficient Descent Advisor (EDA). The second analysis 
was similar but considered the time to fly a specified 
distance to the meter fix. The effects of the test ma- 
trix variables together with the accuracy requirements 
for the predictor will determine the allowable error for 
the predictor inputs. For the B737, the EDA predic- 
tion of meter fix crossing time agreed well with the 
FMS; but its prediction of TOD location probably was 
not sufficiently accurate to enable idle-thrust descents 
in congested airspace, even though the FMS and EDA 
gave si mi lar shapes for TOD location as a function of 
the test matrix variables. For the B777, the FMS and 
EDA gave different shapes for the TOD location func- 
tion, and the EDA prediction of the TOD location is 
not accurate enough to fully enable the concept. Fur- 
thermore, the differences between the FMS and EDA 
predictions of meter fix crossing time for the B777 in- 
dicated that at least one of them was not sufficiently 
accurate. 

Introduction 

In congested airspace today, controllers direct air- 
craft to descend in steps. Since air density, and hence 
drag, increase as the aircraft descends, significant re- 
ductions in fuel consumption and emissions would re- 
sult if aircraft stayed at cruise altitude longer and then 


descended smoothly at idle thrust. The flight manage- 
ment system (FMS) on a large jet can compute the lo- 
cation of top of descent (TOD) assuming an idle-thrust 
descent. To merge aircraft, however, controllers im- 
pose level flight segments, which make it much eas- 
ier for them to estimate the relative speeds of two 
aircraft given their calibrated airspeed (CAS). Three- 
Dimensional Path Arrival Management (3DPAM) en- 
ables arriving aircraft to fly optimized descents com- 
puted by the FMS in congested airspace [1], The 
3DPAM ground automation must accurately predict 
descent trajectories to ensure there are no traffic con- 
flicts and to provide situation awareness to controllers. 
In addition to predicting the meter fix arrival time ac- 
curately to ensure lateral separation, accurate predic- 
tion of the vertical profile is essential to ensure vertical 
separation from aircraft at different altitudes, including 
crossing traffic. 

In the 3DPAM concept, the FMS will be given the 
descent speed profile and the speed and altitude at the 
unknown TOD location and at the meter fix. Assum- 
ing idle thrust in the descent, the FMS will compute 
the vertical profile and fly it with minimal human in- 
tervention. In this study, such descents were flown in a 
simulator with a commercial FMS and no human inter- 
vention. The test matrix varied aircraft weight, descent 
speed, and wind conditions to determine their effect on 
FMS prediction in a B737-700 simulator and a B777- 
200 simulator. 

Due to the potential fuel savings and emissions 
reductions, enabling continuous descents is being pur- 
sued by several research groups [2,3,4,5,6]. Most of 
the previous error analysis has focused on arrival time 
at a waypoint. A few, such as [2,5], have also con- 
sidered the operational TOD prediction error. Using 
a simulator in the current research made it possible to 
use a test matrix, which simplified analyzing and visu- 
alizing the effects of individual factors. Tong, Boyle, 
and Warren [3] and Herndon et al. [7,8] also used sim- 
ulators with commercial FMSs but analyzed the effects 
of different factors; the former considered the de-icing 
setting (which is probably also important for 3DPAM) 
and the latter emphasized comparing aircraft types. 

The first analysis reported in this paper deter- 



mined the effect of the test matrix par ameters on TOD 
location computed by both the FMS and the Efficient 
Descent Advisor (EDA), which is the proposed pre- 
dictor for 3DPAM. The polynomial approximation ap- 
proach was very similar to a previous paper by the 
same author [9]. Analysis of actual TOD locations 
extracted from operational data [10] indicated that the 
TOD path distance can be approximated well by a lin- 
eal - function in cruise and descent altitudes, descent 
CAS, wind, and aircraft weight. This prompted a sec- 
ond look at the FMS predictions of TOD location, fo- 
cusing this time on models linear in descent speed and 
aircraft weight. To better investigate the applicability 
of a lineal - model to multiple aircraft types, the B777 
data previously omitted were included this time. Fur- 
thermore, this paper includes results for predictions 
from EDA, whereas the previous paper only analyzed 
the FMS-computed TOD locations. To interpret what 
is significant variation, previous controller feedback 
suggests that acceptable ground automation TOD pre- 
diction error is likely to be no larger than 5 nmi in order 
to enable optimized descents in congested airspace. 

The second analysis in this paper was similar to 
the first but considered the time to fly a specified dis- 
tance to the meter fix. To interpret these results, con- 
sider that the minimum separation requirement at the 
meter fix is 5 nmi, which is roughly 60 sec in zero wind 
at the meter fix speed and altitude constraints. Conse- 
quently, the ground automation probably needs to pre- 
dict the meter fix arrival time within about 20 sec. 

Since the FMS-computed TOD location should be 
within a few miles of its actual location, comparison of 
the EDA and FMS predictions of TOD location should 
give a reasonable estimate of the operational accuracy 
of EDA, but the data have the advantage of being much 
less noisy than operational data. The FMS prediction 
of the meter fix crossing time, on the other hand, may 
be significantly different from the actual time due to 
failure to maintain the target descent speed or due to er- 
ror in the wind conditions entered into the FMS. Com- 
parison of the EDA and FMS predictions of these times 
can still be useful, however, in identifying significant 
differences between the predictors. Furthermore, the 
accuracy of both predictions will be important in the 
future when data link is used. From both analyses, the 
effects of the test matrix variables will help determine 
the allowable error for the predictor inputs. 

Experiment Procedures 

The experiments in this study were run in sim- 
ulators — a B737-700 and a B777-200 — operated 


by Boeing Phantom Works. These simulators were 
custom-built and each included a commercial FMS, 
which had different manufacturers for these two air- 
craft types. These simulators were designed for de- 
veloping procedures, not pilot training. They include 
random noise. The FMS used both lateral and vertical 
navigation. 

To compute its descent trajectory, the FMS used 
the constraints shown in Figure 1 , which uniquely de- 
termined the trajectory. The parameters shown in each 
box were specified for that segment or that point of the 
descent. In brief, the FMS assumed idle thrust from 
TOD to the meter fix, and it used the speed profile, 
which will be specified in 3DPAM by the controller for 
both cruise and descent. The FMS also used two meter 
fix constraints: altitude 1 1 ,000 ft and speed 250 kt CAS 
(KCAS). In Figure 1, the first segment in the descent 
had constant Mach number, which was the same as the 
cruise Mach number. As the altitude decreased, the 
CAS increased until it reached the target descent CAS, 
which will be specified by the controller in 3DPAM. 
The next segment was then flown at that CAS. The fi- 
nal segment only existed if the descent CAS was faster 
than the meter fix speed constraint, in which case the 
aircraft pitched up to decelerate to that constraint. The 
preceding assumed descent CAS was faster than the 
cruise CAS. If this was not the case, the aircraft first 
decelerated at idle thrust while maintaining level flight 
at cruise altitude or close to it, and then it pitched down 
to maintain the target descent speed at idle thrust. 

The horizontal trajectory was the same nearly- 
straight standard arrival route in all simulation runs. 
Each run started at least 100 nmi from the meter fix 
in cruise at 37,000 ft. To set up the descent in the 
FMS in these experiments, the waypoints were en- 
tered, “ECON MODE” was selected to specify an idle- 
thrust descent, and the FMS-computed descent speed 
was overridden by manually entering the desired de- 
scent speed, which will be specified by the controller 
in 3DPAM. The cruise Mach number was determined 
by the FMS based on a user-specified parameter, but 
it was generally between 0.79 and 0.80. The simula- 
tion runs continued below the meter fix, but analysis 
included only data to the meter fix. 

Two different sets of conditions were used, shown 
in Table 2 and Table 3. The former describes nine runs, 
while the latter describes 22 runs. Five of the entries in 
the two sets were duplicates and run twice. The values 
of the aircraft weight depended upon aircraft type as 
shown in Table 1 . 
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Figure 1. Idle-thrust Descent Schematic 


Table 1. Test Matrix Aircraft Weights 



B737 

B777 

LIGHT_WT 

92,000 lb 

312,0001b 

N OMINAL_WT 

117,000 lb 

379,000 lb 

HEAVY_WT 

131,200 lb 

447,000 lb 


Table 2. First Set of Test Conditions 

Zero Winds in All Nine Runs 


250 KCAS 
LIGHT_WT 

250 KCAS 
N OMINAL_WT 

250 KCAS 
HEAVY_WT 

280 KCAS 
LIGHT_WT 

280 KCAS 
N OMINAL_WT 

280 KCAS 
HEAVY_WT 

320 KCAS 
LIGHT_WT 

320 KCAS 
N OMINAL_WT 

320 KCAS 
HEAVY_WT 


All the trajectories are described by the schematic 
in Figure 1 except those with descent speed 250 KCAS, 
which is slower than the cruise speed. The test condi- 
tions were chosen to cover most of the expected oper- 
ating conditions for descent speed and aircraft weight. 
The distribution of winds aloft and their forecast er- 
rors in the operational environment were unknown, so 
they were chosen somewhat arbitrarily. Since the hori- 
zontal trajectory is nearly straight, the wind fields used 
were actually constant in space and time with direction 
nearly parallel to the arrival route. The wind condi- 
tions labeled “in FMS” means the data entered into the 
descent wind page of the FMS. 


Data 

The analysis of FMS predictions primarily used 
Automatic Dependent Surveillance-Contract (ADS-C) 
messages that Boeing provided. The time between 
ADS-C messages generally did not exceed 64 sec, but 
sometimes it was almost 200 sec. The only ADS-C 
data used in this analysis were the current position and 
the FMS trajectory predictions, which arc in the part 
of the message labeled Intermediate Projected Intent 
Group (IPIG). Each IPIG contains the distance and true 
track from the previous group or from the current po- 
sition for the first group, the predicted barometric alti- 
tude, and the predicted time to reach the point from the 
current position. The path distance of each prediction 
point is computed for analysis by projecting the ini- 
tial position onto the arrival route, computing its path 
distance relative to the meter fix, and adding the given 
distances between waypoints. 

In analysis of the FMS-computed TOD location, 
only the last message sent before the aircraft descended 
from cruise altitude was used. The IPIGs in this mes- 
sage include both the FMS-computed TOD and the me- 
ter fix. In analysis of the FMS-computed time Tj to 
fly distance Df to the meter fix, only the last message 
sent at least distance D j before the meter fix was used. 
Since this prediction rarely starts exactly distance D j 
from the meter fix, linear interpolation was used to es- 
timate Tf. 




Table 3. Second Set of Test Conditions 


Subset A 


Subset B 


Subset C 


117,0001b, 
Zero Winds, 
and One 
Speed Below 


280 KCAS, 
Zero Winds, 
and One 
Weight Below 


280 KCAS, 
117,0001b, 
and One Set of 
Wind Conditions Below 


250 KCAS 
260 KCAS 
270 KCAS 
280 KCAS 
290 KCAS 
300 KCAS 
310 KCAS 
320 KCAS 

Zero Winds in Simulator, 25 kt Tailwind in FMS 
Zero Winds in Simulator, 50 kt Tailwind in FMS 
Zero Winds in Simulator, 25 kt Headwind in FMS 
Zero Winds in Simulator, 50 kt Headwind in FMS 


92,000 lb 25 kt Tailwind in Simulator and FMS 

131 ,200 lb 50 kt Tailwind in Simulator and FMS 

25 kt Headwind in Simulator and FMS 
50 kt Headwind in Simulator and FMS 
25 kt Tailwind in Simulator, Zero Winds in FMS 
50 kt Tailwind in Simulator, Zero Winds in FMS 
25 kt Headwind in Simulator, Zero Winds in FMS 
50 kt Headwind in Simulator, Zero Winds in FMS 


When creating the EDA predictions used in this 
analysis, the appropriate aircraft weight in Table 1 was 
used, even though EDA typically uses the same weight 
for each aircraft type. The cruise Mach number is es- 
timated by the mean of the recorded Mach number in 
cruise from the simulator state data. The initial condi- 
tion for each prediction is at a point with path distance 
211 nmi from the meter fix. For the descent rate for 
the final segment before the meter fix in Figure 1 , the 
data did not indicate an obvious constraint in terms of 
vertical speed, descent rate in ft/nmi, or flight path an- 
gle. Therefore, the EDA predictor simply decelerates 
at level flight, which moves TOD closer to the meter 
fix than it would be if the segment were not level — 
but most likely by less than 1 nmi. 

Analysis of B737 Data 

The rest of this paper discusses the results of the 
data analysis. To make the story easier to follow, this 
section presents all the results for the B737, and the 
next section presents all the results for the B777. The 
details of the analysis arc presented in this section 
along with the results, whereas the section on the B777 
focuses on the results. 


TOD Location 

In this section, the separate effects of descent 
speed, weight, and wind on TOD location arc investi- 
gated. The FMS prediction of the TOD location is crit- 
ical because, as long as the FMS flies the descent, the 
actual TOD will be close to the FMS prediction. There- 
fore, even if the ground automation trajectory predictor 
has more accurate inputs or models than does the FMS, 
the ground prediction of TOD will rarely be closer to 
the actual TOD than is the FMS prediction. Further- 
more, the FMS predictions of TOD location in these 
experiments are the same as they would have been in 
an operational environment. To analyze the TOD loca- 
tion, the trajectory is viewed as two-dimensional, with 
horizontal location specified as distance along the path 
relative to the meter fix, which means the TOD loca- 
tion is then a scalar. 

Effect of Descent Speed and Aircraft Weight 

To quantify the effect of descent speed and weight 
on TOD, the TOD path distance for the runs with zero 
wind (all the runs except those in Subset C of Ta- 
ble 3) was approximated by least squares fit to a linear 
polynomial. The markers in Figure 2 show the FMS- 
computed values of TOD path distance, and the lines 
show the values of this approximation for the weights 
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Figure 2. Speed, Weight Effect on B737 FMS TOD 


and descent speeds in these runs. Table 4 shows the co- 
efficients and goodness-of-fit statistics of this approxi- 
mation. The coefficients also indicate the following: 

• For each 10,000 lb increase in weight, TOD is 
about 4.6 nmi farther from the meter fix. 

• Increasing the descent speed by 10 KCAS 
moves TOD about 3.2 nmi closer to the meter 
fix. 

Assuming the operational ranges of weight and descent 
speed would be close to their ranges in the test ma- 


trix, using a nominal value for either weight or descent 
speed would result in errors exceeding 5 nmi even if 
the predictor were otherwise perfect. 

The next question is whether the TOD locations 
predicted by EDA exhibit the same behavior. The plots 
analogous to Figure 2 look essentially the same, but 
the EDA TOD tends to be about 4 nmi closer than the 
FMS value to the meter fix. This is reflected in the 
differences in the coefficients in Table 4. These results 
confirm that the TOD location determined by the equa- 


Table 4. Speed, Weight Coefficients in TOD Linear Approximation 


Term 

B737 

B777 

FMS 

TS 

FMS 

TS 

Constant 

-110 

-110 

-120 

-170 

Descent CAS (kt) 

0.32 

0.30 

0.29 

0.39 

Weight (lbs) 

-0.00046 

-0.00039 

-0.00013 

-0.00012 

R 2 value 

0.988 

0.990 

0.990 

0.969 

RMS residual 

1.1 nmi 

0.92 nmi 

0.92 nmi 

2.0 nmi 

max abs residual 

1.9 nmi 

1.6 nmi 

1.8 nmi 

5.6 nmi 


tions of motion can be essentially linear in descent 
speed and aircraft weight over their ranges in the test 
matrix used here. 

Effect of Wind 

To characterize the effect of wind on TOD loca- 
tion, first consider the six simulator runs in Subset C 
of Table 3 with the same winds used by the FMS and 
the simulator. The circles in Figure 3 show TOD path 
distance for these runs, and the line shows that this dis- 
tance is linear in the wind, with a 25-kt increase in the 
wind speed moving the TOD location about 4.7 nmi 
farther from the meter fix. 



Figure 3. Effect of Wind on B737 FMS TOD 


There is a simple explanation for the relationship 
between TOD and wind. Since the wind in each run is 
constant over space and time and the speed constraints 


are given as airspeed, the parameters in the equations 
of motion relative to the air mass as well as in their 
constraints are the same in these six runs. Therefore, 
the flight path angle relative to the air mass is the same 
function of altitude in these runs, and it is also the same 
function of time relative to the meter fix time. The 
change in TOD path distance between any two of these 
runs is consequently equal to the difference in the dis- 
tance the air - mass moves relative to ground in the time 
required to descend from TOD to the meter fix. Due 
to small differences in the cruise speeds between runs, 
these descent times vary over 690 to 700 sec in these 
runs. Flence, a change of 25 kt in wind speed should 
move TOD roughly (695 sec x 25 kt = 4.8 nmi), which 
agrees well with the regression line in Figure 3. 

Now consider the runs in Subset C of Table 3 in 
which the simulated winds were different than those 
entered into the FMS descent wind page. These are 
indicated by the markers other than the circles in Fig- 
ure 3. The FMS TOD prediction depends upon the 
error in the winds entered into its descent wind page 
because the FMS blends that wind field with the cur- 
rent wind measurement, which in these experiments is 
a simulated pitot tube measurement. The blending al- 
gorithm is proprietary, so we cannot explain how the 
TOD location depends upon wind error. 

For comparison, the EDA predictions use only 
the wind forecast entered into the FMS, not the sim- 
ulated winds. The TOD locations predicted by EDA 
fit a line very well. A 25-kt increase in the wind 
speed moves the TOD location about 4.5 nmi, which is 
slightly smaller than for the FMS-computed value due 
to shorter descent times — 654 to 658 sec — because 
the TOD is closer to the meter fix. 


Meter Fix Crossing Time 

To compare the descent times, it is insufficient to 
consider only the time between TOD and meter fix, 
since TOD location is not the same in all cases. In- 
stead, we want to compare the time to fly some fixed 
distance Df to the meter fix. To investigate the effect 
of the test matrix variables on the entire descent, the 
value of D f should be larger than the farthest TOD 
from the meter fix but close enough to the meter fix 
that the simulator start-up transients have died by that 
point. For the B737, Df is 95 nmi. 

The portion of the trajectory under consideration 
here consists of a cruise segment and a descent seg- 
ment. For a given simulator run, the ground speed 
should be essentially constant — v c , say — over the 
cruise segment since the Mach number, altitude, and 
wind conditions arc constant except for small random 
errors. The ground speed will clearly vary over the de- 
scent segment, but suppose it has mean value v d . Then 
the predicted time Tj to fly distance Df to the meter 
fix can be approximated as 

Df + Stod Stod 

~ 7 k — , (!) 

Vc V d 

where ,5'tod is the TOD path distance relative to the 
meter fix — wich is negative. This expression shall be 


further manipulated in the next two subsections to ex- 
plain the dependence of Tf on descent speed, aircraft 
weight, and wind. 

Effect of Descent Speed and Aircraft Weight 

For the trajectories with zero wind, ,5’tod can be 
approximated well by a linear polynomial in descent 
speed and aircraft weight, and the airspeed is equal to 
ground speed. Furthermore, for descent speed v rang- 
ing over 250 to 300 KCAS, v is essentially linear in 
v with relative error less than 1%. Therefore, if v d is 
approximately proportional to v, then substitution into 
equation (1) and simplification give 

Tf « 70 + jiv + 72m + jsvw + 7 4 i; 2 (2) 

since v c is roughly the same for all the trajectories with 
zero wind. 

Figure 4 and Table 5 essentially confirm this, but 
they omit the interaction term because it had very lit- 
tle effect. They show not only that this is a good ap- 
proximation but also that the second-order term in de- 
scent speed is significant. For the weights and descent 
speeds in these runs, the solid lines in Figure 4 show 
the values of the approximation that is linear in weight 
and quadratic in descent speed. For each 10,000 lb in- 


Table 5. Speed, Weight Coefficients in Fix Time Approximations 



B737 

B777 

Term 

FMS TS 

FMS TS 


linear model 


Constant 
Descent CAS (kt) 
Weight (lbs) 

1300 

-2.0 

0.00061 

1300 

-1.9 

0.00046 

1600 

-2.3 

0.00010 

1900 
-2.8 
3.0 x 10~ 5 

R 2 value 

0.948 

0.949 

0.946 

0.956 

RMS residual 

1 1 sec 

1 1 sec 

14 sec 

15 sec 

max abs residual 

18 sec 

17 sec 

22 sec 

27 sec 


model lineal - in weight and quadratic in descent speed 


Constant 
Descent CAS (kt) 
Weight (lbs) 
CAS 2 

3000 

-13 

0.00061 

0.020 

2900 

-13 

0.00046 

0.020 

3400 

-16 

0.00010 

0.023 

4000 

-18 

3.0 x 10~ 5 
0.026 

R 2 value 

0.995 

0.997 

0.994 

0.996 

RMS residual 

3.4 sec 

2.8 sec 

4.6 sec 

4.3 sec 

max abs residual 

8.1 sec 

7.1 sec 

9.4 sec 

12 sec 
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Figure 4. Speed, Weight Effect on B737 FMS Fix Time 


crease in weight, Tf is about 6.1 sec longer. If the 
values used in these experiments cover the opera- 
tional range of weight, then using a nominal value in 
the middle of this range would cause errors less than 
0.5 x 0.00061 x (131200-92000) = 12 sec. Whether 
this is acceptable would depend upon the error due to 
other sources. 

Again, the analogous plots for the EDA predic- 
tions look essentially the same, but EDA predicts Tf to 
be 1 to 10 sec shorter than the FMS does. 

Effect of Wind 

To characterize the effect of wind on meter fix 
crossing time Tf, first consider the six simulator runs 


with nominal weight and descent speed and with the 
same winds used by the FMS and the simulator. The 
circles in Figure 5 show Tj for these runs, and the 
solid line in this plot shows that this time is approxi- 
mately lineal - in wind speed, with a 25 -kt increase in 
wind speed decreasing Tf by about 46 sec. For these 
six runs, EDA predicts Tf to be 4 to 8 sec shorter than 
the FMS does. 

To see why Tf is roughly linear in along-track 
wind W v , first recall that the time between TOD and 
meter fix is independent of wind. Furthermore, Stod 
is linear in W v , and ground speed v c in cruise is con- 
stant for each descent and equal to V) + W v , where V) 
is roughly the same for all descents. Hence, as in 


wind :error (kt) : 


TOD Location 


FMS along-track wind (kt) 


-25 ; 
0 



Figure 5. Effect of Wind on B737 FMS Fix Time 


equation (1), 


K + 


Df + (•S'tod + oW x 

# + w v 


( 3 ) 


where 5 *|q D is the TOD path distance for zero wind 

and Vc is the cruise ground speed. For the FMS pre- 
dictions for the B737, K is 695 sec from the above dis- 
cussion of the effect of wind on TOD location, ,S'|q D is 
-75 nmi and a is -0.19 hr from Figure 3, and a plot of 
ground speed for descents with zero wind shows that 
Fc 0) is 454 kt. The dashed line in Figure 5 shows this 
approximation, which suggests it is more accurate than 
the least squares line. 

The markers other than the circles in Figure 5 in- 
dicate the runs in Subset C of Table 3 in which the sim- 
ulated winds were different than those entered into the 
FMS descent wind page. As with the TOD location, 
the effect of this error is not quantified, but compari- 
son with Figure 3 shows that its effect on Tf is roughly 
proportional to its effect on TOD location. 


Analysis of B777 Data 

Now the B777 data is analyzed. The steps in the 
analysis are the same as for the B737 above. 


The first step is again to investigate the separate 
effects of descent speed, weight, and wind on TOD lo- 
cation. 


Effect of Descent Speed and Aircraft Weight 

Figure 6 shows the effect of descent speed and air- 
craft weight on the TOD computed by the FMS for the 
runs in which both the FMS and simulator used zero 
winds. The lines show the values for a linear approx- 
imation, and Table 4 above shows the coefficients and 
goodness-of-fit statistics. The coefficients indicate the 
following: 

• For each 10,000 lb increase in weight, TOD is 
about 1.3 nmi farther from the meter fix. 

• Increasing the descent speed by 10 KCAS 
moves TOD about 2.9 nmi closer to the meter 
fix. 

While the coefficient of aircraft weight for the B777 
is significantly less than for the B737, the range of 
weights is larger. Consequently, the difference in TOD 
location between the extrema of weight is about the 
same for the two aircraft types. As with the B737, us- 
ing a nominal value for either weight or descent speed 
would result in errors exceeding 5 nmi even if the pre- 
dictor were otherwise perfect. 

The first issue with the B777 data is that the linear 
approximation of the EDA predictions of TOD location 
is a significantly worse fit than for the FMS-computed 
values for the B777 as shown by the solid lines in Fig- 
ure 7. Of course, the EDA prediction error for the B777 
is much worse than for the B737 as shown in Figure 8, 
which includes all runs except those in Subset C of Ta- 
ble 3. Not only is the absolute value of the median 
error larger but the variance is also larger. It is not 
known why TOD location as a function of weight and 
descent speed has a significantly different shape in this 
case than it has in the three other cases, but clearly the 
EDA predictor cannot be used for sensitivity studies in 
place of a B777 FMS. 
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Figure 6. Speed, Weight Effect on B777 FMS TOD 
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Figure 7. Speed, Weight Effect on B777 EDA TOD 



Figure 8. Difference Between FMS and EDA TOD 


Effect of Wind 

Figure 9 is the analogue for the B777 of Figure 3. 
For the six simulator runs in Subset C of Table 3 with 
the same winds used by the FMS and the simulator, 
which are represented by the circles in Figure 9, the 
TOD path distance is again linear in wind speed, with 
a 25-kt increase in wind speed moving the TOD loca- 
tion about 5.8 nmi farther from the meter fix. 



Figure 9. Effect of Wind on B777 FMS TOD 


Now consider the runs in Subset C of Table 3 in 
which the simulated winds were different than those 
entered into the FMS descent wind page. These are 
indicated by the markers other than the circles in Fig- 
ure 9. The error in the winds entered into the FMS 


descent wind page seems to have a larger effect than 
for the B737. Furthermore, this wind error is strongly 
correlated with the FMS-computed time between TOD 
and the meter fix, which varies over the range 780 to 
86 1 sec when it should be constant as explained above 
for the B737. This raises questions about how the B777 
FMS is predicting time to fly the trajectory. 

Meter Fix Crossing Time 

Initial transients in the B 777 simulator winds had 
a significant impact on the FMS predictions. For the 
runs in Table 3 , Df cannot be larger than 105 nmi when 
analyzing the FMS estimates. Using this value, Fig- 
ure 10 and Table 5 show that an approximation linear in 
weight and quadratic in descent speed gives a good fit 
of Tf. Again, adding the interaction term has little ef- 
fect. For each 10,000 lb increase in weight, Tf is about 
1.0 sec longer. If the values used in these experiments 
covers the operational range of weight, then using a 
nominal value in the middle of this range would cause 
errors less than 0.5 x 0.00010 x (447000 — 312000) = 
7 sec. Whether this is acceptable would depend upon 
the error due to other sources. 

The EDA predictions of TOD location for the 
B777 are significantly farther from the meter fix than 
the FMS-computed TOD locations, so Df needs to be 
115 nmi when analyzing the EDA predictions. The 
analogous plot to Figure 10 for the EDA predictions 
looks similar, although the variation due to aircraft 
weight is noticeably smaller. Values of Tf based on 
different path lengths D j cannot be compared. Fortu- 
nately, Df can be increased to 1 15 nmi for the FMS es- 
timates for the runs in Table 2 because they started far- 
ther from the meter fix and the initial transients in the 
simulator winds had ended by this point. The resulting 
differences in Tf are strongly correlated with descent 
speed. EDA gives Tf about 35 sec longer than does 
the FMS for descent speed 250 KCAS, 15 sec longer 
for 280 KCAS, and 5 sec shorter for 320 KCAS. The 
reason is unknown. 

Effect of Wind 

Figure 1 1 is the analogue for the B777 of Figure 5. 
For the six simulator runs with nominal weight and de- 
scent speed and with the same winds used by the FMS 
and the simulator, which are represented by the circles 
in Figure 1 1 , Tf is again roughly linear in wind speed. 
The least squares fit line shown by the solid line indi- 
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Figure 10. Speed, Weight Effect on B777 FMS Fix Time 



Figure 11. Effect of Wind on B777 FMS Fix Time 


cates that a 25 -kt increase in wind speed decreasing Tj 
by about 50 sec. The dashed line shows the accuracy 
of equation (3), where for these descents K is 820 sec 
from the above discussion of the effect of wind on TOD 
location, Sj ,q D is -90 nmi and a is -0.23 hr from Fig- 
ure 9, and a plot of ground speed for descents 
with zero wind and nominal weight shows that vi° ! is 
467 kt. Comparison with EDA predictions for these six 
runs is not possible due to the necessary difference in 

Df- 

The markers other than the circles in Figure 1 1 
indicate the runs in Subset C of Table 3 in which the 
simulated winds were different than those entered into 
the FMS descent wind page. As with the TOD loca- 
tion, the effect of this error is not quantified, but com- 
parison with Figure 9 indicates that its relative effect 
on Tf is somewhat less than its relative effect on TOD 
location. 

Conclusions 

The B737 TOD path distance is roughly linear in 
descent speed, aircraft weight, and wind speed. The 
meter fix crossing time is roughly linear in aircraft 
weight and wind speed. The slopes of these lines and 
the accuracy requirement for the predictor will deter- 
mine the allowable error in the inputs used by EDA. 
To meet the current target of error in TOD location 
less than 5 nmi, EDA will need accurate values of de- 
scent speed, aircraft weight, and wind speed. To meet 
the current target of error in meter fix crossing time 
less than 20 sec, EDA will need accurate values of de- 
scent speed and wind speed but weight might not be 


required. The results concerning the effect of wind can 
be generalized to non-constant wind fields [9]. 

Even with perfect inputs, the EDA predictions of 
TOD location for the B737 are 2 to 6 nmi from the 
FMS-computed locations. The cause of this will likely 
need to be identified and reduced in order to enable op- 
timized descents in congested airspace. The fact that 
the EDA and FMS predictions of meter fix crossing 
time agree well for the B737 increases confidence that 
both are using similar equations. 

For the B777, the FMS-computed TOD path dis- 
tance is also roughly linear in descent speed and air- 
craft weight, and using a nominal value for either of 
these would result in absolute error in TOD location 
exceeding 5 nmi. The error in the EDA predictions 
of TOD location is 7 to 18 nmi, and the EDA pre- 
diction of TOD path distance is not linear in descent 
speed and aircraft weight. As with the B737, the most 
likely cause of this is differences between the EDA and 
FMS models of thrust and drag. Furthermore, the dif- 
ferences between the EDA and FMS estimates of me- 
ter fix crossing time indicate that at least one of them 
is not sufficiently accurate. EDA might be the one in 
error, especially given the unexpected shape of TOD 
location as a function of descent speed and weight. On 
the other hand, the B777 FMS estimates of time are un- 
der suspicion since the FMS-computed times between 
TOD and the meter fix depend upon wind conditions 
when theoretical analysis of the equations of motion 
show they should not. 

These experiments only used two aircraft types. 
Stell [10] analyzed operational data for the Airbus 320 
family and the B757. To keep absolute error in TOD 
location less than 5 nmi, the ground predictor must use 
accurate values of descent speed and weight for all four 
of these aircraft types, although nominal weight might 
be acceptable for the B757. Even with these inputs, 
the EDA predictions of TOD location are probably not 
sufficiently accurate for any of these aircraft types to 
fully enable the 3DPAM concept. Therefore, either 
the ground predictor must be improved or the concept 
must be modified. Either approach will require consid- 
erably more data — both laboratory and operational — 
to quantify the prediction errors and validate the con- 
cept. 
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